AD-A233  939 


^89  ”26fl43 

© 

NASA  Technical  Memorandum  101619 
AVSCOM  Technical  Memorandum  89-B-006 


/ 


Scaling  Effects  in  the  Static  Large 
Defection  Response  of  Graphite- 
Epoxy  Composite  Beams 


Karen  E.  Jackson 
Edwin  L.  Fasanella 


L-.s 


fj'ISi 


June  1989 


N/\S/\ 

National  Aeronautics  and 
Space  Administration 

Langley  Research  Center 

Hampton,  Virginia  23665-5225 


SYSTEMS  COMMAND 


AVIATION  RAT  ACTIVITY 


22  0T8 


91  3 


SCALING  EFFECTS  IN  THE  STATIC  LARGE  DEFLECTION  RESPONSE 
OF  GRAPHITE -EPOXY  COMPOSITE  BEAMS 


Karen  E.  Jackson 
Research  Scientist 

U.S.  Army  Aerostructures  Directorate 
NASA  Langley  Research  Center 
Hampton,  Virginia 


Edwin  L.  Fasanella 
Supervisor,  Structures  Support  Unit 
Planning  Research  Corporation 
NASA  Langley  Research  Center 
Hampton,  Virginia 


ABSTRACT 


i 


Scaling  effects  in  the  large 
deflection  response  of  graphite -epoxy 
composite  beams  was  investigated. 
Eight  different  scale  model  beams 
ranging  from  1/6  to  full  scale  were 
subjected  to  an  eccentric  axial  com¬ 
pressive  load  to  promote  large  bending 
deformations  and  failures.  Beams 
having  laminate  stacking  sequences 
including  unidirectional,  angle  ply, 
cross  ply,  and  quasi-isotropic  were 
tested  to  examine  a  wide  "  iety  of 
composite  response  and  fail  modes. 
The  model  beams  were  loaded  under 
scaled  test  conditions  until 
catastrophic  failure.  Data  acquired 
included  load,  end  displacement,  and 
strain  measurements,  and  qualitative 
failure  measurements.  The  experimen¬ 
tal  data  is  compared  to  a  large 
rotation  beam  analysis  and  a  finite 
element  model  analysis  Results  from 
the  tests  indicate  tp'at  the  beam 
response  scales  in  the  small  deflec¬ 
tion  region,  but  deviates  as  the 
response  becomes  nonlinear.  Failure 
modes  are  consistent  between  scale 
models  within  a  laminate  family, 


however,  a  significant  scale  effect  is 
observed  in  strength  of  the  scaled 
beams.  Small  scale  beams  fail  at 
higher  normalized  load  and  sig¬ 
nificantly  higher  normalized  end 
displacement  levels  than  their  full 
scale  prototypes.  It  is  important 
that  this  phenomenon  be  understood 
before  strength  testing  of  scale  model 
composite  structures  can  be  utilized. 


INTRODUCTION 

Scale  model  technology  represents 
one  method  of  investigating  the  struc¬ 
tural  crashworthiness  of  advanced, 
weight  efficient  composite  aircraft 
components  such  as  beams,  frames,  and 
rings.  Impact  tests  on  replica  models 
of  composite  structures  can  provide  a 
cost  effective  alternative  to  full- 
scale  crash  tests.  In  addition,  scale 
model  tests  can  be  conducted  to  verify 
analysis  techniques,  particularly 
finite  element  analyses.  It  is  impor¬ 
tant,  however,  to  understand  the 
limitations  of  scale  modeling  so  that 
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tests  on  sub-scale  models  will  gener¬ 
ate  valid  data.  Scaling  effects  in 

the  response  and  failure  of  composite 
structures  must  be  characterized 

before  the  technique  can  be  used  to 
full  advantage.  A  series  of  tests 
were  conducted  by  Morton  [1]  to  ex¬ 
amine  scaling  effects  in  the  dynamic 
response  of  transversely  impacted 
composite  beams.  Results  from  those 

tests  indicated  that  classical  scaling 
laws  apply  for  elastic  dynamic 
response,  but  a  size  effect  was  ob¬ 
served  as  the  beams  became  damaged 
under  greater  impact  loads. 

The  objective  of  the  current 
research  is  to  investigate  scaling 
effects  in  the  static  large  deflection 
response  of  composite  beams.  The 
scaled  beams  are  loaded  in  a  beam- 
column  fashion  ly  an  eccentric  axial 
compressive  load.  This  testing  con¬ 
figuration  produces  large  bending 
deformations  and  promotes  global 
failure  of  the  beams  away  from  the 
supported  ends.  A  dimensional 

analysis  was  performed  on  the  beam- 
column  system  using  methods  outlined 
in  Baker  [2]  to  determine  the  non- 
dimensional  parameters  or  Pi  terms 
which  govern  the  scaled  response.  An 
experimental  program  designed  to 
validate  the  scaling  laws  was  per¬ 
formed  and  initial  results  are 
reported  in  this  paper.  Also,  a  one 
dimensional  large  rotation  analysis 
and  a  DYnamic.  Crash  Analysis  of 
STructures  (DYCAST)  [3]  finite  element 
model  of  the  composite  beam  were 
developed  for  comparison  with  ex¬ 
perimental  results.  The  results 
obtained  from  the  static  experiments 
'il'iiig  with  a  verified  DYCAST  model 
i  I  I  be  used  to  develop  a  test  matrix 
for  conducting  impact  tests  and 
dynamic  analyses  in  the  future. 

EXPERIMENTAL  PROGRAM 

Beams  having  unidirectional,  angle 
ply,  cross  ply,  and  quasi  -  isotropic 
laminate  stacking  sequences  were 
constructed  of  a  high  modulus  graphite 
fiber  and  an  epoxy  matrix  system 


designated  as  AS4/3502  for  the  static 
tests.  The  full  scale  beam  was  3 
inches  wide  with  a  30  inch  gage  length 
and  48  plies  thick  with  an  average  ply 
thickness  of  0.0054  inches.  The  scale 
model  beams  were  constructed  by  apply¬ 
ing  seven  different  geometric  scale 
factors  including  1/6,  1/4,  1/3,  1/2, 
2/3,  3/4,  and  5/6,  to  the  full  scale 
beam  dimensions.  A  set  of  scaled 
beams  is  illustrated  in  Figure  1  and 
the  dimensions  and  lay-ups  of  each 
beam  are  listed  in  Tabic  1.  The 
thickness  dimension  was  scaled  by 
reducing  the  number  of  layers  in  each 
angular  ply  group  of  the  full  scale 
laminate  stacking  sequence  which 
consisted  of  at  least  six  plies  of 
similar  orientation.  Using  this 
approach,  it  was  not  possible  to 
fabricate  a  1/4  or  3/4  scale  quasi¬ 
isotropic  beam.  Three  replicate 
tests  were  conducted  for  each  laminate 
type  and  size  of  beam.  The  beams  were 
machined  from  panels  which  were  con¬ 
structed  by  hand  from  pre-preg  tape 
and  cured  according  to  manufacturer's 
specifications.  Slight  variations 
were  observed  in  the  thickness  dimen¬ 
sions  of  the  cured  beam  specimens. 
Generally,  the  1/6  scale  beam  was 
thicker  on  a  per  ply  basis  than  the 
full  scale  beam  for  all  laminate 
types .  The  maximum  deviation  in 
normalized  thickness  was  approximately 
six  per  cent . 

During  the  tests  each  beam 
specimen  was  gripped  in  a  set  of 
hinges  which  offset  the  axial  load 
with  a  moderate  eccentricity,  as  shown 
in  Figure  2.  Eight  sets  of  hinges 
were  constructed  to  ensure  that  the 


Identification  of  commercial  products 
and  companies  in  this  paper  is  used  to 
describe  adequately  the  test 
materials.  The  identification  of 
these  commercial  products  does  not 
constitute  endorsement,  expressed  or 
implied,  of  such  products  by  the  U.S. 
Army,  the  National  Aeronautics  and 
Space  Administration,  or  the  pub¬ 
lishers  of  these  conference 
proceedings . 
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end  condition  was  properly  scaled  for 
each  test.  The  hinges  were  pinned  to 
the  platens  of  a  standard  load  test 
machine  which  applied  the  compressive 
vertical  load.  The  h i nged - pinned 
connection  allowed  the  beam  to  undergo 
large  rotations  during  deformation. 
Beam  specimens  were  loaded  until 
catastrophic  failure,  defined  as  loss 
of  load  carrying  capability. 

Each  beam  was  instrumented  with 
back-to-back  strain  gages  located  at 
distances  one-quarter  and  two-thirds 
along  the  length  and  with  strain  gage 
rosettes  at  the  midpoint.  Vertical 
load  was  measured  by  a  load  cell 
located  at  the  base  of  the  bottom 
hinge.  End  displacement  was  measured 
by  an  extensiometer  attached  to  the 
platens  of  the  load  test  machine. 
Vertical  load,  end  displacement,  and 
strain  data  were  recorded  using  a 
personal  computer  based  data  acquisi¬ 
tion  system.  The  analog  signals  were 
amplified  and  filtered  prior  to  being 
digitized  and  converted  to  engineering 
units.  Only  the  load  versus  end 
displacement  data  will  be  presented  in 
this  paper. 


ANALYSIS 

A  one  dimensional  large  rotation 
"elastica"  type  solution  was  developed 
to  predict  the  response  of  the  com¬ 
posite  beam-column  under  eccentric 
axial  load.  The  governing  equation 
for  the  beam  was  derived  from  equi¬ 
librium  of  the  forces  and  moments  on  a 
beam  element.  The  exact  expression 
relating  moment  and  curvature  was 
incorporated  in  the  analysis,  thus 
allowing  the  solution  to  predict  large 
rotation  response.  The  solution  of 
the  governing  equation  is  outlined  in 
Timoshenko  and  Gere  [4]  for  the 
"elastica"  problem  and  was  adapted  for 
this  problem  by  applying  the  end 
moment  boundary  conditions  produced  by 
the  eccentric  vertical  load.  The 
solution  is  given  in  terms  of  elliptic 
integrals  and  predicts  the  end  dis¬ 
placement,  transverse  displacement  of 
the  midpoint  of  the  beam,  and  end 


rotation  for  increasing  load. 

The  beam  bending  stiffness  was 
derived  from  the  method  described  by 
Whitney  [5]  in  which  the  bending 
stiffness,  El,  from  classical  beam 
theory  is  replaced  by  an  equivalent 
stiffness  for  the  composite  beam.  The 
beam  is  considered  as  a  special  case 
of  a  laminated  plate  in  which  the 
length  is  much  larger  than  the  width. 
Consequently,  the  transverse  displace¬ 
ment  is  assumed  to  be  a  function  of 
the  axial  coordinate  only.  Also,  only 
the  moment  along  the  axial  direction 
is  assumed  to  be  present.  This  is 
analogous  to  a  plane  stress  assumption 
in  elasticity.  The  transverse  and 
twist  curvatures  are  expressed  in 
terms  of  the  axial  curvature  and  the 
bending  stiffnesses,  and  are  then 
substituted  into  the  equation  relating 
bending  moment  and  axial  curvature . 
The  equivalent  beam  bending  stiffness 
is  the  coefficient  of  the  axial  curva¬ 
ture  from  this  equation.  It 
incorporates  the  shear  and  twist 
coupling  terms  which  are  important  for 
angle  ply  and  quasi- isotropic 
laminates . 

In  addition  to  the  beam  analysis, 
the  nonlinear  finite  element  struc¬ 
tural  analysis  computer  program  DYCAST 
[3]  was  used  to  model  the  composite 
beam- column.  Since  the  DYCAST  program 
will  be  used  to  model  the  beam-column 
under  impact  conditions  in  the  future, 
the  static  case  was  developed  to 
verify  the  model  and  to  compare  with 
the  large  rotation  beam  solution.  The 
composite  laminate  was  discretized 
into  60  beam  elements  which  were 
constrained  to  permit  only  planar 
deformations ,  as  shown  in  Figure  2 . 
The  hinges  at  the  top  and  bottom  of 
the  beam  were  modeled  by  two  rigid 
beam  elements  each.  The  model  assumed 
pinned  conditions  between  the  load 
machine  and  the  hinge,  and  clamped 
conditions  between  the  hinge  and  beam. 
The  bending  stiffness  used  in  the 
DYCAST  model  was  the  same  as  used  in 
the  beam  analysis  outlined  previously. 
The  complete  model  had  192  degrees  of 
freedom.  The  applied  load  was  in¬ 
creased  incrementally  at  one  end  using 
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a  static  full  Newton  iterative  tech¬ 
nique  in  which  the  stiffness  matrix 
was  updated  in  each  iteration.  The 
full  Newton  procedure  was  required 
since  the  modified  Newton  method  which 
updated  the  stiffness  matrix  for  each 
load  step  failed  to  converge  in  the 
nonlinear  region  of  the  response 
curve . 


RESULTS 

Normalized  load  versus  end  dis¬ 
placement  plots  and  corresponding 
photographs  of  a  complete  (1/6  through 
full  scale)  set  of  failed  beam 
specimens  for  the  unidirectional, 
angle  ply  cross  ply,  and  quasi¬ 
isotropic  laminates  are  shown  in 
Figures  3-6.  Vertical  load  was  nor¬ 
malized  by  the  Euler  column  buckling 
load  for  the  beam,  and  end  displace¬ 
ment  was  normalized  by  the  gage 
length.  Since  three  repeat  tests  were 
performed  for  each  laminate  type  and 
size  of  beam,  the  results  from  one 
representative  test  are  presented 
here.  Repeatability  between  the  three 
tests  was  good. 

Normalized _ Load _ Versus _ End 

Displacement  Results 

In  general,  the  load  versus 
displacement  curves  show  that  the 
response  scales  for  small  end  dis¬ 
placement  ratios,  typically  less  than 
0.1.  Deviation  from  scaled  response 
is  observed  for  all  laminate  types  as 
the  beams  undergo  large  deflections 
and  the  response  becomes  nonlinear. 
The  angle  ply  beams  show  the  most 
pronounced  deviation  from  scaled 
response,  as  seen  in  Figure  4(a).  The 
small  scale  beams  fail  at  a  higher 
normalized  load  and  end  displacement 
l-'/ul  than  the  full  scale  beam.  This 
observed  scale  effect  in  failure 
behavior  is  significant.  The  1/6 
scale  beams  fail  at  an  end  displace¬ 
ment.  to  length  ratio  from  2  to  10 
times  the  value  for  the  full  scale 
beam  depending  on  the  laminate  type. 


Failure  Mechanisms 

The  photographs  shown  ir.  Figures 
3(b)  through  6(b)  indicate  that  while 
the  failure  modes  for  the  laminate 
types  considered  in  this  study  are 
different  from  each  other,  they  are 
similar  between  scaled  beams  within 
the  laminate  family.  Failure  modes 
appear  to  be  independent  of  specimen 
size.  The  unidirectional  beams,  shown 
in  Figure  3(b),  failed  by  fiber  frac¬ 
tures  near  the  midpoint  of  the  beam. 
This  failure  mode  is  typical  of  all 
the  unidirectional  beams  r/6  ’"hrough 
full  scale.  Failure  of  the  angle  ply 
beams  occurred  by  transverse  matrix 
cracking  along  45  degree  fiber  lines. 
There  was  no  evidence  of  fiber 
breakage,  as  shown  in  Figure  4(b). 
The  cross  ply  laminates  exhibited 
combined  failure  mechanisms  of 
transverse  matrix  cracking  and  fiber 
fracture.  As  the  cross  ply  beam 
underwent  large  rotations,  the  90 
degree  plies  located  in  the  center  of 
the  laminate  developed  transverse 
matrix  cracks .  The  cracks  were  evenly 
spaced  and  resulted  in  uniform  pieces 
of  debris,  some  of  which  are  shown  in 
Figure  5(b)  for  the  5/6  scale  beam. 
The  ultimate  failure  of  the  cross  ply 
beam  was  caused  by  fiber  fractures  in 
the  0  degree  plies.  The  quasi¬ 
isotropic  beams  failed  through  a 
combination  of  matrix  cracking, 
delamination,  and  some  fiber  failure. 
Although  the  photograph  in  Figure  6(b) 
does  not  give  a  good  indication,  the 
damaged  quasi - isotropic  beams  are 
highly  curved.  The  sequence  of 
failure  events  occurred  such  that  the 
remaining  intact  section  of  the  beam 
consisted  of  an  unsymmetric  laminate, 
resulting  in  the  observed  curvature. 

Analytical  Results 

Comparison  of  the  experimental 
data  for  the  1/6  and  full  scale 
specimens  with  the  large  rotation  beam 
analysis  and  the  DYCAST  finite  element 
analysis  is  plotted  in  Figures  7(a) 
through  7(d)  for  each  of  the  laminate 
types.  Agreement  between  the  two 
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analysis  methods  is  excellent,  even 
though  they  approach  the  problem  in 
different  manners.  The  beam  solution 
assumes  an  inextensible  beam,  while 
the  DYCAST  model  allows  in- plane 
deformations  due  to  membrane  loads. 
Also,  the  large  rotation  beam  analysis 
incorporates  the  exact  nonlinear 
expression  for  beam  curvature,  while 
the  DYCAST  model  uses  the  linear 
expression.  However,  these  factors  do 
not  appear  to  be  important  in  the 
response  prediction.  DYCAST  appears 
to  be  sufficiently  accurate  for  future 
dynamic  analyses  of  the  scaled  beam 
since  a  closed  form  "elastica"  type 
analysis  is  not  available. 

Good  correlation  is  obtained 
between  the  experiment  and  the  beam 
solution  and  DYCAST  for  small  load 
ratios,  generally  less  than  0.4. 
However,  both  the  large  rotation  beam 
solution  and  the  DYCAST  model  typi¬ 
cally  overpredict  the  experimental 
beam  response  as  the  load  ratio  and 
normalized  end  displacement  values 
increase  and  beam  rotations  become 
large.  This  is  true  for  all  of  the 
laminates  tested.  The  slope  of  the 
response  curve  in  the  large  deflection 
region  (normalized  end  displacement 
greater  than  0.2)  as  predicted  by  both 
analyses  is  in  good  agreement  with 
experiment,  as  shown  in  Figures  7(a), 
7(c),  and  7(d).  Overprediction  of  the 
response  by  the  large  rotation  beam 
analysis  and  DYCAST  may  be  due  in  part 
fo  certain  assumptions  made  in  the 
analysis  including  constant  stiffness 
assumptions.  The  stiffness  of  the 
beam  is  reduced  due  to  nonlinear 
matei ial  properties  and  damage  events 
■  ai  It  as  l  i  aiisvi- rse  in. it  i  i  x  crocking, 
which  are  not  modeled  by  the  ana  lysis. 

DISCUSSION 

The  results  presented  here  indi¬ 
cate  that  a  significant  scale  effect 
exists  in  the  failure  behavior  since 
the  smaller  scale  beams  fail  at  a 
higher  normalized  load  and  much  higher 
normalized  end  displacement  value  than 
the  full  scale  beam.  Stress  and 


strain  based  failure  criterion  such  as 
maximum  stress,  maximum  strain,  Tsai- 
Hill,  or  Tsai-Wu,  would  not  be  able  to 
predict  the  observed  scale  effect. 
According  to  classical  scaling  laws, 
stress  and  strain  should  scale  as 
unity.  Consequently,  under  perfectly 
scaled  experimental  test  conditions 
the  stress  and  strain  in  a  model  beam 
will  be  the  same  as  for  the  prototype. 
A  stress  analysis  of  the  scaled  test 
will  predict  one  value  of  end  dis¬ 
placement  to  length  ratio  at  which 

failure  should  occur.  Morton  [1] 
discusses  a  linear  elastic  fracture 

mechanics  approach  to  the  strength 
scaling  of  transversely  impacted 
composite  beams  and  shows  that  a 
theory  for  a  notch- sensitive  or 
brittle  material  can  predict  scaling 
effects  in  a  cracked  plate. 

Application  of  these  theories  to  a 
stress  analysis  of  the  beam-column 
problem  is  planned  as  a  continuation 
of  the  experimental  and  analytical 
results  presented  here. 

CONCLUDING  REMARKS 

Scaling  effects  in  the  large 
deflection  response  and  failure  be¬ 
havior  of  graphite -epoxy  composite 
beams  was  investigated.  A  series  of 
static  tests  on  scale  model  composite 
beams  having  unidirectional,  angle 
ply,  cross  ply,  and  quasi-isotropic 
laminate  stacking  sequences  was 
conducted.  The  beams  were  loaded 
under  an  eccentric  axial  compressive 
load  to  promote  large  bending  deforma¬ 
tions  and  global  failure.  Plots  of 
norma  I  i zed  load  versus  end  displace 
inent  wore  generated  to  compare  with  a 
one  dimensional  large  rotation  com¬ 
posite  beam  analysis  and  a  DYCAST 
finite  element  model. 

Results  from  the  experiments  show 
that  beam  response  scales  in  the  small 
deflection,  elastic  region;  however, 
deviations  from  scaled  response  appear 
as  the  beams  undergo  large  deflections 
and  rotations.  The  degree  of  varia¬ 
tion  from  scaled  response  is  dependent 
on  laminate  stacking  sequence.  Angle 
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Table  1.  Scale  model  beam  test  specimen  dimensions  and  lay-ups . 
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Figure  3.  Unidirectional  graphite-epoxy  composite  beam  results 
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(a)  Normalized  load  versus  end  displacement. 
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(a)  Normalized  load  versus  end  displacement. 
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(b)  Failed  beam  specimens. 

Figure  6.  Quasi-isotropic  graphite-epoxy  composite  beam  results. 
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(a)  Unidirectional  beam  results. 
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(b)  Angle  ply  beam  results. 

Figure  7.  Comparison  of  normalized  load  versus  end  displacement 
experimental  data  for  1/6  and  full  scale  beams  with  DYCAST  finite  element 
analysis  and  large  rotation  exact  solution. 
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(c)  Cross  ply  beam  results. 


(d)  Quasi-isotrcpic  beam  results. 

Figure  7.  Comparison  of  normalized  load  versus  end  displacement 
experimental  data  for  1/6  and  full  scale  beams  with  DYCAST  finite  element 
analysis  and  large  rotation  exact  solution. 
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16  Abstract 


Scaling  effects  in  the  large  deflection  response  of  graphite- 
epoxy  composite  beams  was  investigated.  Eight  different  scale  model 
beams  ranging  from  1/6  to  full-scale  were  subjected  to  an  eccentric 
axial  compressive  load  to  promote  large  bending  deformations  and 
failures.  Beams  having  laminate  stacking  sequences  including  uni¬ 
directional,  angle  ply,  cross  ply,  and  quasi-isotropic  were  tested 
to  examine  a  wide  variety  of  composite  response  and  failure  modes. 

The  model  beams  were  loaded  under  scaled  test  conditions  until  cat¬ 
astrophic  failure.  Data  acquired  included  load,  end  displacement, 
and  strain  measurements,  and  qualitative  failure  measurements.  The 
experimental  data  is  compared  to  a  large  rotation  beam  analysis  and 
a  tinite  element  model  analysis.  Results  from  the  tests  indicate 
that  the  beam  response  becomes  nonlinear.  Failure  modes  are  con¬ 
sistent  between  scale  models  within  a  laminate  family,  however,  a 
significant  scale  effect  is  observed  in  strength  of  the  scaled  beams. 
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